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The misalignment between the orbital plane of a transiting exoplanet and the spin 
axis of its host star provides important insights into the system's dynamical history. 
The amplitude and asymmetry of the radial-velocity distortion during a planetary 
transit (the Rossiter-McLaughlin effect) depend on the projected stellar rotation rate 
C"| vsml and misalignment angle A, where the stellar rotation axis is inclined at angle 

Oh I to the line of sight. The parameters derived from modelling the R-M effect have, 

q however, been found to be prone to systematic errors arising from the time-variable 

5^ asymmetry of the stellar spectral lines during transit. Here we present a direct method 

for isolating the component of the starlight blocked by a planet as it transits the 
host star, and apply it to spectra of the bright transiting planet HD 189733b. We 
1-1 model the global shape of the stellar cross-correlation function as the convolution of 

a limb-darkened rotation profile and a gaussian representing the Doppler core of the 
average photospheric line profile. The light blocked by the planet during the transit is 
I a gaussian of the same intrinsic width, whose trajectory across the line profile yields 

\q a precise measure of the misalignment angle and an independent measure of vsml. 

We show that even when vsml is less than the width of the intrinsic line profile, 
If} the travelling Doppler "shadow" cast by the planet creates an identifiable distortion 

in the line profiles which is amenable to direct modelling. Direct measurement of 
the trajectory of the missing starlight yields self-consistent measures of the projected 
stellar rotation rate, the intrinsic width of the mean local photospheric line profile, 
the projected spin-orbit misalignment angle, and the system's centre-of-mass velocity, 
r* Combined with the photometric rotation period, the results give a geometrical measure 

. of the stellar radius which agrees closely with values obtained from high-precision 

transit photometry if a small amount of differential rotation is present in the stellar 
photosphere. 
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1 INTRODUCTION 



The recent discoveries of strong projected spin-orbit mis- 
alignments in the transiting exoplanets XO-3b (jHebrard et| 
al.|2008||Winn et al.|2009a) , HD 8 0606b (|Winn et al.|2009bT 



Pont et al.|2009[ ) and WASP-14b ( |Johnson et al.|2009|, and 



retrograde orbital motion in WASP-17b ( jAnderson et al.| 
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0828(A). The data are publicly available at the CDS (CDSARC.U- 
strasbg.fr). 
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2009[ | and HAT-P-7b ( |Winn et al.|2009c||Narita et al.|2009 1 

indicate violent dynamical histories for a significant fraction 
of the population of hot- Jupiter planets. 

These misalignments are measured by obtaining 
densely-sampled radial-velocity observations during tran- 
sits. The starlight blocked by the planet during a transit 
possesses the radial velocity of the obscured part of the stel- 
lar surface. The radial- velocity centroid of the light emanat- 
ing from the visible parts of the stellar disk exhibits a reflex 
motion disturbance -the Rossiter-McLaughlin (RM) effect 
- whose form depends on the projected stellar equatorial 
velocity vsini, the impact parameter b of the planet's path 
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across the stellar disc, the projected misalignment angle A 
between the orbital axis and the stellar spin axis, and the 
relative sizes of the star and planet. 

Detailed semi-analytic formulations for determining the 
anomalous departure of the line centroid from the stellar re- 
flex orbit have been develo ped by [Ohta, Taruya, fc Suto| 
(20051), iGimenezI (120061) andlHirano et al.l d2009b. The pro- 



cess of measuring and interpreting this reflex velocity shift is 
not, however, straightforward. If the stellar rotation profile 
is even partially resolved by the spectrograph, the missing 
light within the planet's silhouette introduces an asymme- 
try into the line spread function. For instruments such as 
SOPHIE and HARPS, radial velocities are derived using a 
gaussian fit to a cross-correlation function (CCF) computed 
from the stellar spectrum and a weighted line mask ( |Baranne| 
et al.|1996 Pepe et al.|2002 \ . Such methods work extremely 
well for the symmetric and time-invariant profile shapes en- 
countered outside transit. When applied to profiles with a 
time-varying degree of intrinsic asymmetry, however, they 
yield velocity measures that depart systematically from the 
velocity centroid of the visible starlight. 



Winn et al. ( 2005 1 overcame this difficulty for their 



iodine-cell observations of HD 209458b using a direct mod- 
elling approach which they have employed in subsequent pa- 
pers. They built models of both the unobscured stellar pro- 
file and the light blocked by the planet into their analysis 
of the line-spread function, deriving semi-empirical correc- 
tions to the model velocities to enable meaningful compari- 
son with the data. Nonetheless, their study of the RM effect 
in HD 189733 ( |Winn et aL]|2006[ ) shows a clear pattern of 
correlated radial-velocity residuals during the transit. In a 
more recent HARPS study of the RM effect during a transit 



of HD 189733b, Triaud et al. ( 2009 1 noted an almost iden- 



tical systematic pattern of correlated residuals between the 
radial velocities of the gaussian fits to the CCFs, and the 
line centroid velocities of the best-fitting model computed 



using the formulation of Gimenez ( 2006 1 



In developing an analytic model of this velocity anomaly 



for cross-correlation spectra, Hirano et al. ( 2009 1 have shown 



that such errors are exacerbated as v sin I increases and 
the asymmetry becomes more extreme. If uncorrected, the 
anomaly can lead to significant over-estimation of the value 
of vsinl, particularly for rapidly-rotating planet-host stars. 
The motivation for the present study is to develop a com- 
prehensive model of the changes in CCF morphology that 
occur before, during and after a transit, in terms of the pro- 
jected stellar equatorial rotation speed, the spin-orbit mis- 
alignment angle, the intrinsic width of the combined stellar 
CCF and instrumental broadening function, and the stel- 
lar limb-darkening coefficient. In principle this approach is 



similar to that of Winn et al. (20051, but we use it to de- 



compose the observed CCF into its various components and 
track them directly. A similar methodology has already been 
used to model the RM effect in the binary stars V1143 Cyg 
( jAlbrecht et al.|2007[ ) and DI Her |Albrecht et al.|p)09"| ). To 
illustrate the effectiveness of this approach as an alternative 
to gaussian-fitting in the study of transiting exoplanets, we 
re-analyse the same HARPS observations of HD 189733b 



described by Triaud et al. (2009). We show that the spec- 



tral signature of the light blocked by the planet is clearly 
discernible in the residuals when a model of the unobscured 
starlight is subtracted from the data, and track it directly 



to derive new measurements of the projected stellar rotation 
velocity and spin-orbit misalignment angle. 



2 OBSERVATIONS AND DATA ANALYSIS 

Densely-sampled spectral time-series of HD 189733 were ob- 
tained with the HARPS instrument on the ESO 3.6-m tele- 
scope at La Silla on the nights of 2006 September 8 and 2007 
August 29. The observations are described in detail by |Tri-| 
aud et al. (2009). The 2006 September data were taken with 



a cadence of 10.5 minutes, while the 2007 data were sam- 
pled every 5.5 minutes. The signal-to-noise ratio (SNR) in 
the continuum of each spectrum was typically 230 in echelle 
orders 61 through 70 of all spectra in the 2006 sequence, and 
130 for the 2007 data which employed shorter exposures. 

The spectra were extracted and numerical cross- 
correlation functions (CCF) were computed using the 
most recent version of the HARPS data-reduction software 
(Mayor et al. 20091. This involves order- by-order cross- 



correlation of the extracted stellar spectrum with a weighted 



mask function, as described by Baranne et al. ( 1996 ) and 



Pepe et al. (20021. The mask used for these observations 



was derived from the line list for a star of spectral type 
K5. The SNR in the pseudo-continuum to either side of the 
cross-correlation peak was measured to be ~ 6300 for a typ- 
ical spectrum in the 2006 series, and ~ 3900 in the 2007 
sequence. The cross-correlation process thus yields a multi- 
plex gain factor of approximately 29 in SNR, for both data 
sets, allowing the SNR in any individual CCF to be esti- 
mated directly from the SNR of the original spectrum. 



2.1 Stellar CCF profile model 

The CCF of a slowly-rotating star exhibits a flat continuum 
and a single dip, which can be fitted to high precision with 
a gaussian profile. We assume that the local line profile at 
any point on the surface of a limb-darkened, rotating star 
takes the form of a gaussian, 



9{x) 



2ns 



exp f— x 2 /2s 2 ) 



(1) 



where the dimensionless velocity x — v/v sin / and gaussian 
sigma s = a /vsinl are expressed in units of the projected 
stellar equatorial rotation speed v sin I. 

We assume a standard linear limb darkening model with 
limb-darkening coefficient u: 



B(v) = B{1){1 -M + u/i) 



(2) 



using the standard convention that the direction cosine of 
the local surface normal relative to the line of sight is 
H = cos 9, with 6 = at the centre of the stellar disc. Assum- 
ing solid-body rotation, the limb-darkened rotation profile is 
obtained by substituting fi = \J\ — x 2 — y 2 and integrating 
Eq. [2] between limits y = ± y/l — x 2 , where x = v/v sin / as 
before. The resulting expression is normalised by integrating 
over the range — 1 < x < 1 to obtain ( |Gray|1976 1 



fix) 



6 ((1-u) vT 



l)/4) 



7T (3- 
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These expressions are normalised such that 
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Figure 1. Coordinate system used in definition of missing-light 
model. The instantaneous orbital coordinates of the planet, x p 
and z p , are projected on the plane of the sky. The coordinate 
Up denotes projected distance from the stellar rotation axis. In 
this example, the projected obliquity A is a positive angle. At 
this obliquity, the planet's trajectory, with impact parameter b, 
carries it mainly across redshifted parts of the stellar surface, as 
denoted by the colour scale. 



g(x)dx 



f(x)dx ■ 



(4) 



The observed rotation profile is the convolution of the two 
functions, 



K x ) = / f(z)g(x - z)dz, 



(5) 



which is computed via a straightforward numerical integra- 
tion. Since the CCFs are computed in the velocity frame of 
the solar-system barycentre, the shift of the model CCF to 
match a Keplerian orbit solution is achieved by defining Xij 
for a pixel i with velocity Vij in the barycentric frame to be 



(^(ecoso; + cos(uj + w)) + 7) 



(6) 



where Vj is the true anomaly at the time of the jth obser- 
vation and all other orbital elements K, e,ui, 7 have their 
usual meanings. 

The model of the jth observation of the stellar CCF is 
thus h(xij). 
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Figure 2. Residual maps of the time-series CCFs for the night 
of 2006 September 7/8 (upper) and 2007 August 28/29 (lower). 
The model stellar profile h{j has been subtracted, leaving the 
bright travelling signature of the starlight blocked by the planet 
during the transit. This is seen superimposed on a fixed-pattern 
residual spectrum ai which remains stable during each night but 
changes subtly from one night to the next. The horizontal dotted 
line gives the phase of mid transit. Crosses denote times of first, 
second, third and fourth contact. The stellar radial velocity is 
plotted as a dotted curve flanked by dashed lines at ±u sin I. 



2.2 Missing-light model during transit 

At any moment, the position of the planet on the plane of 
the sky is given by 



r sin(^ + to — it/ 2) 

r cos(f + u — 7t/2) cos i 



(7) 
(8) 



where r is the instantaneous distance of the planet from the 
star, and i is the inclination of the orbital angular momen- 
tum vector to the line of sight. The z axis is parallel to the 
projected direction of the orbital axis (Fig. [T]). If x p , z p and 
r are expressed in units of the stellar radius, the fraction of 



the starlight blocked by the planet during the total part of 
a transit is 



R p 1 — u + u/j, 
M l-u/3 



(9) 



where fi = y/ 1 — x\ — z| . 

The radial velocity of the missing starlight relative to 
the star's instantaneous velocity depends on the projected 
distance of the planet from the stellar rotation axis. Follow- 
ing |Winn et al.| ( |2005| l we define A = <^ sp in — or t>it where <fi 
denotes position angle in the plane of the sky. 
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Figure 3. Residual maps of the time-series CCFs for the night 
of 2006 September 7/8 (upper) and 2007 August 28/29 (lower). 
The scaled best-fit model hij of the stellar profile, the travelling 
planet signature j3 g(xij — u p ) and the fixed sidelobe pattern ai 
have all been subtracted from the data as in Eq. |16| leaving only 
noise. 



Rotating the coordinate frame in the plane of the sky, 
we obtain the perpendicular distance of the planet from the 
stellar rotation axis in units of the stellar radius: 

Up = x p cos A — z p sin A, (10) 

which in turn yields the barycentric velocity of the missing 
light, v p = Up v sin/. 

The combined model of the stellar CCF and the missing 
starlight during the total part of the transit is thus 



Pij = h(xij) + (3 g(xij - Up). 



(11) 



Analytic expressions for computing f3 and u p during the par- 
tial phases of the transit are given by |Ohta, Taruya, fc Suto] 
( 20051) and iGimenezI 12006b. 



The term /3 g(xij—u p ) represents the travelling gaussian 
"bump" produced by the blocked starlight. The appearance 



of this residual bump after subtraction of the model stellar 
profile from the observed CCFs is illustrated in greyscale 
form in Fig. [2] The coefficient j3 depends on the ratio R p /R* 
and the limb- darkening coefficient u. The width of the gaus- 
sian g gives an additional constraint on the non-rotating 
local CCF width s, given the reasonable assumption that 
Rp « R*. The trajectory of the planet signature through 
the line profile is described by u p , which in turn depends on 
wsin/ and A. These separate constraints on s, u and vsinl 
are valuable, since they help to break any degeneracy in 
the description of the shape of the stellar CCF, particularly 
when vsinl is small compared to the intrinsic linewidth. 



2.3 Model fitting and error scaling 

The data and model are both orthogonalised by subtracting 
their optimal mean values: 



d' 



= Pa 



(12) 
(13) 



using inverse- variance weights uiij — 1/oy. The multiplica- 
tive constant A by which the model must be scaled to give 
an optimal fit to the data is then given by the expression 



(14) 



The CCF obtained by cross-correlating a late-type stel- 
lar spectrum with a line mask inevitably contains weak side- 
lobes, caused by random alignments between stellar lines 
and mask lines at arbitrary shifts. The cross-correlation dip 
therefore appears superimposed on a background pseudo- 
continuum which exhibits weak fixed-pattern residuals. 
These appear as an alternating pattern of light and dark 
vertical stripes in Fig. [2] Imperfect modelling and removal 
of the spectrograph blaze function prior to cross-correlation 
can also produce a weak asymmetry in the CCF profile, 
which would leave a similar pattern in the residuals. At this 
stage in the calculation, the sidelobe pattern is removed by 
computing and subtracting the optimal average of the resid- 
ual spectra: 



(15) 



The sidelobe pattern can also be influenced by moonlight 
and telluric absorption levels, and potentially by weak, time- 
variable asymmetries in the CCF caused by starspots. This 
produces subtle changes in the sidelobe pattern from one 
night to the next. Little change is seen within a single night. 
The method is therefore most effective when the spectrum- 
averaging summations are performed separately for different 
nights. The residuals following the subtraction of the side- 
lobes and the signature of the light blocked by the planet 
are shown in Fig. [3] 

The goodness of the model fit to the observed CCF is 
then measured via the statistic 



= ^2(d'ij - Ap'ij - a l fw l] . 



(16) 
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Table 1. System parameters from recent studies of HD 189733b. Values adopted from previous work arc italicised. 
Parameter 



Winn et al. 



(2006) 



Pont et al. 



(2007) 



Triaud et al. 



I 2009 I 



This study 



Units 



P 

K s 
e 



i 

M, 
R, 

Rp/ i?* 



2453937.7759 ± 0.0001 2453931.12048 ± 0.00002 2453988.80339 ± 0.00006 2453988.8034 



2.218575 




86.1 ± 0.2 
0.82 ± 0.03 
0.73 ±0.02 



2.218581 ± 0.000002 


0.671 ± 0.008 
85.68 ± 0.04 
0.82 ± 0.03 
0.75 ±0.01 
0.1572 ±0.0004 



2.218573 ± 0.000001 
0.202 ±0.001 
0.0041 ± 0.0022 
-24 ± 34 
0.687 ± .006 
85.51 ± 0.1 
0.82 ± 0.03 
0.77 ± 0.01 
0.1581 ± 0.0005 



2.218573 
0.202 
0.0041 
-24 
0.687 
85.51 
0.82 
0.77 
0.1581 



d 

km s _1 
degrees 



M© 

Re 



Care is needed when estimating the error bars on both the 
data and the fitted parameters. The spectral resolving power 
of HARPS is R = A/AA = 120000, yielding a velocity res- 
olution of 2.5 km s _1 . The CCFs produced by the HARPS 
data reduction pipeline are oversampled, being binned in 
velocity increments of 0.25 km s _1 . The errors on adjacent 
points in the CCF are therefore strongly correlated. In order 
to ensure that all data points contributing to \ 2 are statis- 
tically independent, we bin the data and the fitted model in 
velocity, by a factor ten before computing \ 2 ■ Initially we 
estimate the relative uncertainties ay of the CCF data in 
each spectrum from the SNR of the original spectrum using 
the multiplex gain factor described in Section [2] above. 



2.4 Markov-chain Monte Carlo parameter fitting 

Given a sequence of CCFs derived from a set of spectra 
densely sampled during transit and distributed around the 
rest of the orbit, the overall \ 2 statistic is therefore a func- 
tion of the orbital parameters K, e, ui and 7; the orbital 
period P and the epoch T p of periastron; the parameters 
v sin 7, s and u describing the shape of the rotationally- 
broadened stellar CCF and the continuum tilt; and the error 
scaling factor /. 

The limited phase range of the HARPS observations 
does not permit as reliable a determination of the full 
set of orbital parameters as the combined SOPHIE and 
HARPS data used by |Triaud et~aL 



(2009). We therefore 



chose not to fit the full radial-velocity orbit and the photo- 
metric transit parameters simultaneously with the spectro- 
scopic transit parameters. Instead, we fixed the values of all 
other model parameters at the values given by |Triaud et al.| 
(2009) in their recent combined photometric and spectro- 



scopic analysis of the same system. These values are sum- 
marised and compared with their values from earlier studies 
in Table [l] Only the system centre-of-mass velocity, 7, was 
re-determined directly from the data. This was done because 
the star is known to exhibit starspot activity, so small night- 
to-night changes in 7 are to be expected as the star rotates. 
Although spot activity could in principle affect the appar- 
ent value of K too, the latter is well-determined from RV 
observations collected over many orbits. Any degradation of 
the fit by starpot-induced changes in the apparent radial 
acceleration during the transit is expected to be small. 

We determined the model parameters vsml, A, s and 
7 and their joint posterior probability distribution using a 
Markov-chain Monte Carlo (MCMC) approach. 



Strongly-correlated pairs of parameters are best avoided 
in MCMC, since they lead to slow convergence and unde- 
sirably strong correlations between successive steps in the 



chain (Tegmark et al. 2004). We expect the parameters s 



and v sin I, which together determine the width of the stellar 
CCF outside transit, to be correlated to some extent, since 
their quadrature sum is fixed by the width of the observed 
CCF. To maintain orthogonality between jump parameters 
we replaced s with the full width at half maximum (FWHM) 
of the CCF, which is approximated by 



cf = [(vsml)' 2 ±«g] 1/2 = vsin/[l + 4s 2 In 2] 



1/2 



(17) 



Here v g is the FWHM of the gaussian representing the local 
stellar and instrumental profile: 



2s vs'ml^ 



(18) 



The Markov chain is constructed by repeatedly comput- 
ing x 2 for a sequence of parameter sets {v sin/, A, v cc f , 7}. 
At each step k, a set of trial parameters is computed by 
applying gaussian perturbations to the preceding set: 

usin/ fc = «sin/ fc _i + fa vsin iG(0, 1) 
A fc = Afc_i + fa\G(0, 1) 

Vccf,k = «ccf,fe-l + /0"„ ccf G(0, 1) 

7fc = 7fe-i+MG(0,l) 

The scale factor / is of order unity. The goodness-of-fit 
is recalculated for the new parameter set, and the change 
Ax 2 = xt ~ xt-i i s computed. The new set is either ac- 
cepted or rejected according to the Metropolis-Hastings rule: 
if Ax 2 < the set is accepted, and if Ax 2 > 0, the new step 
is accepted with probability e _Ax I 2 . If the proposed pa- 
rameters are accepted, they are written to the next record 
in the Markov chain. If not, the previous set of values is 
recorded. The scale factor / is tuned to ensure that roughly 
25 percent of all steps are accepted once the chain has set- 
tled down to a steady state; by trial-and-error we found that 
/ = 0.7 achieved the desired acceptance rate. 

The chain evolved to a steady state within about 100 ac- 
cepted proposals. We terminated the initial "burn-in" phase 
after at least 100 proposals had been accepted, when the 
value of Xfe exceeded the median x 2 °f au previous steps 
for the first time ( |Knutson et al.||2008[ ). At this stage we 
refined our estimates of the variances of the binned CCF 
data, replacing the initial estimates erf,- with <xfjXmin,.j l n d- 
Here Xminj ls the contribution of the jth spectrum to the x 2 
statistic for the best-fitting model in the chain so far, and 
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Figure 4. Correlation diagrams for the joint posterior probability distributions of the four model fitting parameters {tisin/, A, v s , 7} for 
the transit of HD 189733b observed on 2006 September 7/8. 




Figure 5. Correlation diagrams for the joint posterior probability distributions of the four model fitting parameters {tisin/, A, v s , 7} for 
the transits of HD 189733b observed on 2007 August 28/29. 



The Doppler shadow of HD 189733b 7 



rid is the number of degrees of freedom associated with the 
spectrum. This ensures that \ 2 is scaled correctly, and that 
the posterior probability distributions of the fitted parame- 
ters have the correct variances. 

We then ran the chain for a further 100 proposal accep- 
tances, and re-evaluated the variances of the four fitting pa- 
rameters from the chains themselves. After a further short 
burn-in period, we performed a production run of 10 4 ac- 
cepted proposals. The correlation lengths of all four chains 
were computed as described by |Tegmark et al.| |2004 ) , and 
found to be of order 10 steps in all four parameters. 

The profile decomposition method carries higher com- 
putational overheads than the anomalous RV approaches 



developed by Winn et al. ( 2005 1 and Hirano et al. ( 2009 1 



Substantial speed gains can be realised by performing the 
numerical integration to generate the synthetic stellar profile 
only once per step in the Markov chain, and shifting the pro- 
file to the appropriate orbital velocity at the time of each ob- 
servation. The MCMC analyses presented here, with chains 
comprising 10 4 accepted steps, incurred execution times of 
order 50 minutes on a standard 2.66 GHz laptop processor, 
as opposed to a few minutes for a conventional analysis. 



3 RESULTS 

In Figs.|4]and[5]we show the correlation diagrams for all four 
jump parameters, for the nights of 2006 September 7/8 and 
2007 August 28/29. In these plots we show the FWHM u g of 
the intrinsic gaussian profile, rather than the FWHM w cc f of 
the convolved profile. Contrary to our initial expectations, 
we find that v sin / and u g are only very weakly correlated. 
This confirms that the two independent measures of both 
quantities (from the shape of the stellar profile and the width 
and trajectory of the missing-light signature) serve to fix 
their values unambiguously. 

The four fitted parameters for each of the two transits 
are listed in Table [2] for three different values of the linear 
limb-darkening coefficient u. These values bracket the value 



u ~ 0.725 from the ATLAS 9 models of Claret (2000) for a 



stellar effective temperature of 5050K, loggr = 4.5 and near- 
solar metallicity, as reported by |Bouchy et al.| ( |2005[ ) for HD 
189733. The tabulated values illustrate that the best-fitting 
values of vsinl and v s are only very weakly dependent on 
the precise value of u. 

Although the two sets of observations were made nearly 
a year apart, the fitted values of the four parameters agree 
within their uncertainty ranges. The radial velocity 7 of 
the system's centre of mass remained constant to within 
10 m s" 1 . The spin-orbit misalignment angle A is found to 
be slightly negative in both years. We find an absolute value 
slightly, but not significantly, smaller than the (0.85±0.30)° 
measured by Triaud et al. The uncertainty in A is also 0.3 
degrees, comparable with the uncertainty derived by |Triaud| 
( 2009 1 from the same data. 

1 is signifi- 



et al. 



Our value of v sin / 



3.10 ± 0.03 km s" 
-1 



Triaud 



cantly lower than the 3.321q'q7 km s" 1 obtained by 
|et al.| ( |2009| ) using an uncorrected fit to the HARPS ra- 
dial velocities. As Triaud et al. noted, the radial-velocity 
measurements obtained during transit by gaussian fitting 
show systematic departures from the intensity-weighted av- 
erage velocity of the unobscured parts of the stellar disc. 



Table 2. Mean and one-sigma errors for MCMC model parame- 
ters fitted to the transit events recorded on 2006 September 7/8 
and 2007 August 28/29. Results are shown for three MCMC anal- 
yses using different values of the linear limb-darkening coefficient 



Parameter 2006 Sept 7/8 2007 Aug 28/29 



Units 



u = 0.675 
v sin I 
A 

7 

u = 0.725 
v sin I 
A 

"g 
7 

u = 0.775 
v sin I 
A 

7 



3.071 ±0.019 
-0.37 ±0.234 
6.277 ± 0.024 
-2.209 ±0.010 



3.0835 ± 0.020 
-0.40 ±0.24 
6.286 ± 0.022 
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They modelled the departure of the radial- velocity centroid 
from the gaussian-fitting velocity as a function of v sin I, and 
found that vsinl = 3.05 km s _1 yielded the best fit to the 
observed pattern of residuals. An identical pattern of residu- 
als was seen in the earlier R-M study of this system by | Winn] 



et al. (20061. Their empirical correction procedure yielded 
a lower v sin / = 2.97 ± 0.22 km s -1 . Our direct tracking of 
the missing-starlight component in the CCF eliminates the 
need for such empirical corrections to indirect measures of 
the missing-light velocity, and yields a more precise result 
that lies within the error ranges of both studies. 



4 DISCUSSION AND CONCLUSIONS 

In principle we could determine a lower limit on the stel- 
lar radius from v sin I and the photometric rotation pe- 
riod, which was determined by Henry & Winn ( 2008 1 to 



be 11.953 ±0.009 days. The very close alignment of the pro- 
jected orbital and spin axes in the plane of the sky makes 
it highly probable that the inclination I of the stellar spin 
axis to the line of sight should be very close to the incli- 
nation i = 85.5 degrees of the planet's orbit. Assuming 
sin J = sini = 0.997, our value for «sinl yields a direct 
estimate of the stellar radius _R* = 0.732 ± 0.007 R©, in- 
dependently of the usual assumptions involving isochrone 
fits to the stellar density and effective temperature. This is 
significantly smaller than the R, = 0.77 ± 0.01 R0 deter- 



mined by Triaud et al. (20091, and marginally less than the 
R* = 0.75 ±0.01 R Q found by |Pont et al"1 ( |2007 |. 

We suspect that the discrepancy arises from differential 



rotation of the stellar photosphere. Gaudi & Winn (20071 



pointed out that under some circumstances, the RM effect 
could be used to measure stellar differential rotation. As 
with most previous studies of this kind, we have measured 
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v sin I under the assumption that the star rotates as a solid 
body. In all stars for which differential rotation has been 
measured directly, however, a solar-like pattern is found: 
the equator rotates faster than the poles, with the depar- 
ture from the equatorial spin rate having a sine-squared de- 
pendence on latitude. The typical difference in rotational 
angular frequency AS! = f2 G quator — ^poic is found to be of 
order 0.04 radian per day for rapidly-rotating stars K dwarfs 
similar to HD 189733, and only weakly dependent on rota- 
tion rate (Barnes et al. 20051. This is very similar to the 



ened stellar spectrum; the result would then be used as the 
reference spectrum. A very similar procedure was described 



by Winn et al. ( 2005 I for calibrating the departure of the 



velocities measured during a transit of HD 209458b from the 



solar value of Af2. 

If the starspots giving rise to the photometric modu- 
lation signal are concentrated at the same stellar latitude 
as the planet's path across the stellar disc, we should derive 
the correct value of the stellar radius. Unfortunately we have 
no way of determining the dominant latitude from which the 
modulation signal originates. If the orbital and spin axes are 
closely aligned, the impact parameter of the planet's path 
across the star carries it across a stellar latitude of 39° in the 
hemisphere facing away from the observer. Although |Pont| 
et al. (20071 noted the passage of the planet across one or 



two isolated dark spots, these spots are significantly fore- 
shortened. It is reasonable to expect the modulation signal 
to arise mainly from spots at intermediate to low latitude in 
the observer's hemisphere. A difference in latitude of 20 or 
30 degrees is sufficient to give a difference in angular veloc- 
ity 8Q. ~ 0.015 radians per day between the main spot belts 
and the latitudes traversed by the planet. This discrepancy 
yields a stellar radius that is too small by about 3 percent, 
which is sufficient to reconcile our radius estimate with those 



of |Pont et al.| ( [2007| > and |Triaud et al.| |2009). 

We conclude that direct modelling of the cross- 
correlation function during a planetary transit is feasible 
even for a host star such as HD 189733 whose vsml is com- 
parable to the intrinsic line width. The profile decomposi- 
tion method described here yields values and error estimates 
for the stellar spin rate and orbital obliquity that agree 
closely with the results of previous studies of this system. 
The parameter values are free of the systematic errors that 
occur when velocity measurements are derived from gaus- 
sian fits to line profiles with inherent time-variable asym- 
metry. This circumvents the need for semi-empirical correc- 
tions when modelling the Rossiter-McLaughlin effect. The 
method makes the fullest use of all the information present 
in the cross-correlation function: v sin / is tightly constrained 
by the shape of the stellar rotation profile as well as by the 
spectral signature of the starlight blocked by the planet. The 
stellar radius obtained from the photometric rotation period 
and the stellar v sin I agrees well with values obtained from 
high-precision transit photometry, if a modest amount of 
differential rotation is present. 

The method described here could in principle be ex- 
tended to observations for which an iodine cell is used to 
track the spectrograph PSF and wavelength scale. The anal- 
ysis of such observations requires a reference spectrum of 
the host star, taken with the same instrument without the 
iodine cell. An artificial reference spectrum could be gen- 
erated by convolving a limb-darkened rotation profile with 
the the spectrum of a narrow-lined star of the same spec- 
tral type. The Doppler signature of the light blocked by 
the planet could be mimicked by scaling, shifting and sub- 
tracting the same narrow-lined spectrum from the broad- 



predictions of the Ohta, Taruya, & Suto ( 2005 \ model, using 



high-resolution solar spectra. 
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